1 



0^ 



Antiferromagnetism in four dimensions: search for non-triviality* 

J. L. Alonso'^, H. G. Ballesteros^, I. Campos'^, J. M. Carmona'^, J. Clemente Gallardo'', 

L. A. Fernandez'', D. Iniguez^, V. Martm-Mayor'', A. Munoz Sudupe'', A. Tarancon** and C. L. UUod*^ 

'^Departamento de Fi'sica Teorica, Facultad de Ciencias, 
Universidad de Zaragoza, 50009 Zaragoza, Spain 

''Departamento de Fi'sica Teorica I, Facultad de CC. Fi'sicas, 
Universidad Complutense de Madrid, 28040 Madrid, Spain 



o 



We present antiferromagnetism as a mechanism capable of modifying substantially the phase diagram and the 
critical behaviour of statistical mechanical models. This is particularly relevant in four dimensions, due to the 
connection between second order transition points and the continuum limit as a quantum field theory. 

We study three models with an antiferromagnetic interaction: the Ising and the 0(4) Models with a second 
neighbour negative coupling, and the RP^ Model. Different conclusions are obtained depending on the model. 
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An antiferromagnetic (AFM) system shows in 
general different critical behaviour from that of its 
ferromagnetic counterpart, as wiU be seen below. 
This can be used to study the triviality problem in 
Quantum Field Theory, by computing the critical 
exponents or the renormalized coupling, which is 
expressed in terms of the Binder cumulant C/l of 
the model as 
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and should be different from zero in a non-trivial 
theory. 

One way by which one could hope to see new 
critical behaviour, i.e. new universality classes, 
is by considering AFM models. Actually, it has 
been shown in low dimensions that antiferromag- 
netism can produce new critical behaviour. For 
example, in the Ising Model we can introduce 
it by means of a second-neighbour (related with 
higher derivatives in the action) negative cou- 
pling, which produces frustration in the system. 
This makes a new ground state to appear, which 
in 2d gives a new transition line with new critical 
behaviour 

The 3d RP^ AFM Model has also recently been 
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studied M . The Hamiltonian is given by 



(2) 
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where Vi is a three component normalized vec- 
tor. Here a negative coupling produces a mixed 
ground state, formed by one sub-lattice oriented 
in a fixed direction and the other lying on the 
plane perpendicular to it, which brings about a 
new universality class. 

We have considered the following models in 4c?. 

1. AFM Ising Model 

This is the simplest AFM model: the Ising 
Model in a hypercubic lattice with first- and 
second-neighbour couplings 



(3) 



It was simulated on a y = lattice with periodic 
boundary conditions, using a Heat Bath update 
(clusters algorithms are not effective in most of 
AFM models). The Finite Size Scaling (FSS) was 
studied up to L = 24 |]. 

This model turns out to have a rich phase dia- 
gram, with different vacuum structures. There is 
a region of ferromagnetic order (FM): Un = oo 
(fJo stands for a fixed spin); there is also a 
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Figure 1. Phase Diagram of the Ad AFM Ising 
Model. 



phase, HPAF, which consists of a FM config- 
uration on a three-dimensional cube and AFM 
on the other direction /z: o-„ = (— 1)"*'(To; fi- 
nally, in the FAF phase, we have FM on a two- 
dimensional plane and AFM on the other two di- 
rections: an = (— l)"f'+"''(To, where /x, i/ < v) 
can be any directions. We define order parame- 
ters according to these structures. 

The phase diagram is shown in Fig. [I| PM 
means a disordered phase and the dotted lines 
are the transition lines given by Mean Field the- 
ory. The transition line connected with the usual 
Ising point (/?2 = 0) is second order with the same 
(classical) exponents. The PM-HPAF line is first 
order. And we find another line, PM-PAF, which 
is disconnected from the Ising critical point. This 
transition is clearly first order for values of f3i 
greater than 0.2, but the behaviour is not so clear 
for /3i < 0.2. 

We produced a great amount of data at the 
points /3i — 0.1, 0.05, 0. At /?i the system de- 
couples into two independent sublattices, known 
as F4 lattices in the literature. We found that the 
two peaks characteristic of a first order transition 
are only distinguished aX (3i — 0, and with the 
maximum size, L = 24. 

This AFM model teaches us that we have to be 
very careful, because first order transitions with 
a very large correlation length can deceive you if 
you do not go to sufficiently large lattice sizes to 
see the real behaviour. With smaller sizes, appar- 



ent second order transitions arise, and even false 
critical exponents can be measured. This system 
behaves for small L as a weak first order tran- 
sition, and presents the pseudo critical exponent 
a/iy — 1, precursor of a first order behaviour for 
large lattice sizes [Q. 

2. AFM RP2 Model 

This is the four-dimensional version of the 
model defined in (^. It was simulated by com- 
bining Metropolis and Over-relaxed algorithms, 
up to L = 24 lattice sizes 

Two order parameters appear due to the special 
nature of the ground state, one FM and the other 
AFM (or staggered). They are defined through 
tensors attached to every lattice site 

T^f = vfvf-l6-0, (4) 

as the normal or staggered sum of these tensors 
to the whole volume. 

In order to calculate critical exponents, we use 
the FSS ansatz that allows to write 

_ (0(21,13)) _ /, Fo(^(2£,/3)/2£) 

(0(L,f3)) Fo(aL,P)/L) ^ ^ 

(except for scaling corrections), so that 

Qo|q,.2 = 2""/'^ + ---. (6) 

from where the critical exponent xq is readily 
extracted. 

The model has a second order transition at 
(3 « —1.34, its critical exponents being almost 
gaussian, the differences likely due to logarith- 
mic corrections. To know about the triviality of 



Table 1 

5fl(L,/3c(oo)) for the M AFM RP^ Model. 
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~8 3.16 (3) 
10 2.84 (4) 
12 2.61 (5) 
16 2.34 (4) 
20 2.08 (9) 
24 1.95 (13) 
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the theory we should therefore look at the renor- 
malized coupling, defined by (|l[), at the critical 
temperature. As can be seen in Table |, for the 
staggered sector the ga gets smaller values as the 
lattice size grows. We obtain a good fit to 1/ log L 
without a constant term, so that the theory seems 
to be trivial. However, in the ferromagnetic sec- 
tor, gji takes a negative, rather stable, value. This 
might be an indication of a non-trivial limit for 
this theory. 
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Figure 2. Quotients to obtain ^/v (0(4)). 



3. AFM 0(4) Model 

We consider a system of spins {^n} taking val- 
ues in the hyper-sphere S'^ C R"* and placed in the 
nodes of a hypercubic lattice. The interaction is 
defined by the Hamiltonian: 

-(32 ^ ■ 



n 
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The phase structure is found to be qualitatively 
very similar to that of the Ising Model of Fig. ^ 
but now we have a continuum symmetry. We did 
Monte Carlo simulations in lattices ranging from 
L = 6 to L = 24 (1). The update method was 
a combination of Over-relaxed and Heat-Bath al- 
gorithms, being z ~ 1. 

In the 0(4) AFM Model, the whole phase tran- 
sition line PM-PAF seems to be second order. We 
computed the critical exponents on this line at 
/3i = coupling {F4 lattice) with FSS techniques. 
The results are quoted in Table ||. The correla- 
tion length exponent, is very close to the Mean 
Field prediction, namely ly = 0.5. Although we 
have not obtained a reliable measure of the spe- 
cific heat exponent our results point to a/i^ < 1 
and disagree with the expected behaviour for a 
weak first order transition. 



Table 2 

Critical Exponents of the Ad AFM 0(4) Model. 



Lattice sizes 


7/;/ 


(3/,y 


V 


6/12 


2.417(3) 


0.791(4) 


0.474(10) 


8/16 


2.403(3) 


0.792(6) 


0.483(8) 


10/20 


2.410(2) 


0.790(4) 


0.471(6) 


12/24 


2.403(5) 


0.797(5) 


0.483(7) 



From our 7/1^ estimation (see Fig. ^) (or (ijv 
using the hyper-scaling relation) the exponent t] 
associated with the anomalous dimension of the 
field is 77 « —0.4. This fact itself would im- 
ply the non-triviality of the theory, but bearing 
in mind the results of section 1, transitory be- 
haviour cannot be discarded. However, the stabil- 
ity of our measures of 7/1/ for the different lattice 
sizes, (more than a hundred of standard devia- 
tions from 7y = 0), makes this rather unlikely. We 
also measured the Binder cumulant at the crit- 
ical point, finding that it stays almost constant 
when increasing the lattice size, which points to 
a non-zero value of the renormalized constant in 
the thermodynamical limit. 
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